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Periodic Emission from the Gamma-Ray 
Binary 1FGL J1018.6-5856 

The Fermi LAT Collaboration* 

Gamma-ray binaries are stellar systems containing a neutron star or black hole, with gamma-ray 
emission produced by an interaction between the components. These systems are rare, even though 
binary evolution models predict dozens in our Galaxy. A search for gamma-ray binaries with the 
Fermi Large Area Telescope (LAT) shows that 1FGL J1018.6-5856 exhibits intensity and spectral 
modulation with a 16.6-day period. We identified a variable x-ray counterpart, which shows a sharp 
maximum coinciding with maximum gamma-ray emission, as well as an 06V((f)) star optical 
counterpart and a radio counterpart that is also apparently modulated on the orbital period. 1FGL 
] 1018.6-5856 is thus a gamma-ray binary, and its detection suggests the presence of other fainter 
binaries in the Galaxy. 


T wo types of interacting binaries containing 
compact objects are expected to emit gam- 
ma rays (/): microquasars — accreting black 
holes or neutron stars with relativistic jets (2 ) — and 
rotation-powered pulsars interacting with the wind 
of a binary companion (3). Microquasars should 
typically be powerful x-ray sources when active, 
and hence such gamma ray-emitting systems may 
already be known x-ray binaries. Indeed, the bright 
x-ray source Cygnus X-3 is now known to be 
such a source ( 4, 5). The existence of pulsars in- 
teracting with stellar companions of early spectral 
types is predicted as an initial stage in the forma- 
tion of high-mass x-ray binaries (HMXBs) con- 
taining neutron stars (6). These interacting pulsars 
are predicted to be much weaker x-ray emitters 
and may not yet be known or classified x-ray 
sources. Gamma-ray binaries may thus not be as 
rare as they appear to be, and many systems may 
await detection. 

A gamma-ray binary is expected to show 
orbitally modulated gamma-ray emission due 
to a combination of effects, including changes in 
viewing angle and, in eccentric orbits, the degree 
of the binary interaction, both of which depend 
on binary phase. Periodic gamma-ray modulation 
has indeed been seen in LS 5039 (period 3.9 


*All authors with their affiliations appear at the end of this 

paper. 


days), LS I +61° 303 (26.5 days), and Cygnus X-3 
(4.8 hours) (4, 7, 8\ and gamma-ray emission is 
at least orbital phase-dependent for the PSR 
B 1259-63 system (3.4 years) (9). However, the 
putative gamma-ray binary HESS J0632+057, 
for which a 321 -day x-ray period is seen, has not 
yet been shown to exhibit periodic gamma-ray 
emission (10). PSR B 1259-63 contains a pulsar, 
and LS 5039 and LS I +61° 303 are suspected, 
but not proved, to contain pulsars, whereas 
Cygnus X-3 is a black hole candidate. A search 
for periodic modulation of gamma-ray flux from 
LAT sources may thus lead to the detection of 
further gamma-ray binaries, potentially revealing 
the predicted HMXB precursor population. The 
first Fermi LAT (11) catalog of gamma-ray sources 
(“1FGL”) contains 1451 sources (12), a large 
fraction of which do not have confirmed counter- 
parts at other wavelengths and thus are poten- 
tially gamma-ray binaries. 

To search for modulation, we used a weighted 
photon method to generate light curves for all 
1FGL sources in the energy range 0. 1 to 200 GeV 
(13). We then calculated power spectra for all 
sources. From an examination of these, in addi- 
tion to modulation from the known binaries LS 
I +61 ° 303 and LS 5039, we noted the presence 
of a strong signal near a period of 16.6 days from 
1 FGL J1 01 8.6-5856 (Fig. 1). 1FGLJ 10 18.6-5856 
has a cataloged 1- to 100-GeV flux of 2.9 * 10 -8 


photons cm -2 s _l , making it one of the brighter 
LAT sources. The source’s location at right as- 
cension (R.A.) = lO* 1 18.7™, declination (decl.) - 
-58° 56.30' (J2000; ±1.8', 95% uncertainty) means 
that it lies close to the galactic plane (b = -1.7°), 
marking it as a good candidate for a binary sys- 
tem. 1FGL J10 18.6-5856 has been noted to be 
positionally coincident with the supernova rem- 
nant G284.3-1.8 (12) and the TeV source HESS 
J1 018— 589 (14), although it has not been shown 
that these sources are actually related. 

The modulation at a period of 16.6 days has 
a power more than 25 times the mean value of 
the power spectrum and has a false-alarm prob- 
ability of 3 x 1 o~ 8 , taking into account the num- 
ber of statistically independent frequency bins. 
From both the power spectrum itself (15) and 
from fitting the light curve, we derived a period 
of 16.58 ± 0.02 days. The folded light curve 
(Fig. 1) has a sharp peak together with additional 
broader modulation. We modeled this to deter- 
mine the epoch of maximum flux 7^ by fitting 
a function consisting of the sum of a sine wave 
and a Gaussian function, and obtained T max ~ 
modified Julian date (MJD) 55403.3 ± 0.4. 

The gamma-ray spectrum of 1FGL J1018.6- 
5856 shows substantial curvature through the LAT 
passband. To facilitate discussion of the lower- 
energy (<1 GeV) and higher-energy (>1 GeV) 
gamma rays, we adopted as our primary model 
a broken power law with photon indices T 0 . n and 
T j_ 10 for energies below and above 1 GeV, re- 
spectively. The best-fit values (13) are r 0 .i_i = 
2.00 ± 0.04** ± 0.08 and T^q = 3.09 ± 0.06a* ± 
0. 12^*, along with an integral energy flux above 
100 MeVof (2.8 ± 0.1 8tat ± 0.3^0 x 10“ 10 erg 
cm -2 s" 1 . A power law with exponential cutoff 
(7, 8), dN/dE = N 0 (E( GeV)“ r exp (-E/E c ), gives 
an acceptable fit with F = 1.9 ±0.1 and£ c = 2.5± 
0.3 GeV (statistical errors only). Although this 
spectral shape is qualitatively similar to that of 
pulsars and of LS I +61 ° 303 and LS 5039, so far 
no detection of pulsed gamma-ray emission has 
been reported (16). 

To investigate variability on the 16.6-day 
period, we folded the data into 1 0 uniform bins in 
orbital phase and then refit the broken power-law 
parameters within each phase bin. The resulting 


www.sciencemag.org SCIENCE VOL 335 13 JANUARY 2012 


189 


Downloaded from www.sciencemag.org on June 28, 2012 



REPORTS 


folded light curve (Fig. 2) indicates substantial 
variability in both the source brightness and 
spectral shape. In agreement with the detection 
of multiple harmonics of the orbital period in 
the power spectrum, there appear to be two pri- 
mary features. For phases 0.2 to 0.6, the spectral 
curvature decreases and the peak of the spectral 
energy distribution lies below the LAT passband 
(indicated by r 0 .i_i > 2). The onset of this soft 
spectrum is approximately coincident with a rise 
in x-ray emission and a peak in radio emission 
(see below). A weaker peak appears in the low- 
eneigy (< 1 GeV) y-ray flux at phase 0.5 (Fig. 2). 
For the remaining phases, the LAT spectrum 
hardens with a comparatively sharp rise to, and fall 
from, a peak around 1 GeV (T 0 .i_i < 2, T i_i 0 > 2). 
The variable spectral shape implies that only a 
modest fraction of the flux could represent steady 
magnetospheric emission from a pulsar. 

We undertook observations of the location of 
1FGL 11018.6-5856 covering the energy range 
0.3 to 10 keV by means of the X-ray Telescope 
(XRT) onboard the Swift satellite. The first ob- 
servation was obtained on 29 September 2009 
with an exposure of 5 ks. A single source was 
detected in the XRT image (Fig. 3) within the 
LAT error circle. We then obtained additional ob- 
servations torn January to April 2011 to search 
for x-ray variability (73) and found large am- 
plitude variability. Folded on the gamma-ray 
ephemeris (Fig. 4), there is a sharp peak in x-ray 
flux, coincident with the gamma-ray peak. How- 
ever, in addition to this, a sine wave-like periodic 
modulation is also seen that peaks near phase 
0.3 to 0.4. 

Swift Ultraviolet/Optical Telescope (UVOT) 
(17) observations were obtained simultaneously 
with the x-ray observations. The x-ray source is 
positionally coincident with a bright source seen 
in the UVOT images (Fig. 3) (73), which in turn 
is coincident with a source in the U.S. Naval 
Observatory B 1 .0 catalog at (J2000.0) R.A. =* Up 
1 8 m 55 s .60 ± 0.1", decl. =-58° 56' 46.2" ± 0.1 
Spectroscopic observations of the optical coun- 
terpart were performed using the South African 
Astronomical Observatory 1.9-m telescope and 
the 2.5-m telescope at the Las Campanas Obser- 
vatory. Absorption lines due to H, He i, and He u 
identify it as an early-type star. We used a spectral 
atlas (18) to estimate the spectral type. He n 
X4686 is present in absorption, which indicates 
a main-sequence star. The ratio of He n A.4541 to 
He i X4471 implies an 06 spectral type. Weak 
emission is seen from N m but not He n, which 
indicates an ((f)) classification. We therefore es- 
timate the spectral type as 06V((f )). This is very 
similar to the spectral type of LS 5039 (79). 
Interstellar absorption bands provide an estimate 
of the reddening, E(B - V ), defined as the rel- 
ative absorption in the B and V optical bands; 
from the features at 4430 and 5780 A, we de- 
rive E(B -V) = 0.9 and 1 .6, respectively. Taking 
V ~ 12.6 from measurements with the All Sky 
Automated Survey (ASAS) (20), we derive a 
distance, d, to 1FGL J101 8.6-5856 from Earth of 


5 ± 2 kpc, allowing for uncertainties in the 
reddening and spectral classification. 

Radio observations of the 1FGL J1018.6- 
5856 region were obtained with the Australia 
Telescope Compact Array (ATCA) at frequencies 
of 5.5 and 9 GHz. A faint radio source at R.A. = 
10 h 18 m 55 s .580, decl. — 58° 56' 45.5" (±0.1" 
and 0.3", respectively) is coincident with the 
stellar position. The radio source was clearly 
seen to be variable (Fig. 4). Unlike the gamma- 
ray and x-ray modulation, there is no obvious 
brightening in the radio at phase zero. Instead, 
it appears that the radio may be following the 
smoother sine wave-like component of the x-ray 
modulation. 

1FGL J1 01 8.6-5856 shares many properties 
with LS 5039. They are both fairly steady gamma- 
ray sources on long time scales, their periodic 
modulations have not shown large changes, and 
their optical counterparts are of a very similar 
spectral type. The x-ray light curve of LS 5039 
appears to be highly repeatable (21, 22), and the 
x-ray light curve of 1FGL J10 18.6-5856 also 
shows repeatable behavior with a flux increase 
around phase 0 repeated over four orbital periods. 
The lack of variability in UV/optical brightness 
is also reminiscent of LS 5039 (23, 24 ). This 
suggests that there is little ellipsoidal modulation 
of the primary star and hence that it substantially 
underfills its Roche lobe. On the other hand, the 
relative phasing of the gamma-ray spectral mod- 
ulation and flux modulation differ from those of 
LS 5039 where the spectrum is softest when the 


flux is highest (8). Also, for LS 5039 the phases 
of maximum x-ray and gamma-ray do not co- 
incide (8, 22). The brightest peak in the folded 
gamma-ray light curve of 1 FGL J 1 0 1 8.6-5856 at 
phase 0 is associated with the hardest gamma-ray 
spectrum and is coincident with x-ray flaring and 
minimum radio emission. Finally, 1FGL J1018.6— 
5856 has a much longer orbital period. 

The gamma-ray modulation observed in 1FGL 
J101 8.6-5856 could be due to anisotropic in- 
verse Compton (IC) scattering between stellar 
photons and high-energy electrons that varies 
with orbital phase, as proposed for LS 5039 and 
LS I +61° 303 (7, 8). However, the modulation 
amplitude is considerably lower in 1 FGL J1 01 8.6- 
5856 [(/n^-ZminVamax +/mta) * 25%] relative 
to that of LS 5039 (« 60%). Modulation ampli- 
tude should increase with eccentricity and is high- 
est for systems viewed edge-on (25); however, in 
the case of LS I +61° 303, the modulation frac- 
tion has been observed to undergo large changes 
(26). If the IC scattering interpretation is correct, 
then this implies that 1FGL J1 01 8.6-5856 has 
both low inclination and low eccentricity. For com- 
parison, the eccentricity of LS 5039 has been re- 
ported to be in the range of 0.3 to 0.5 (79, 27, 23). 
Although a low inclination angle implies that it 
would be difficult to measure die radial velocity 
of the companion from optical studies, the small 
Doppler shifts predicted would facilitate a pulsa- 
tion search at GeV energies. 

The gamma-ray spectral variability of 1FGL 
J10 1 8.6-5856 over the orbit is also reminiscent of 
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Fig. 1. Power spectrum of the LAT weighted photon light curve (E > 100 MeV) of 1FGL ]1018.6-5856. 
The power spectrum is oversampled by a factor of 4 relative to its nominal resolution. The red dashed line 
indicates the 16.6-day period; the blue dashed lines are the second, third, and fourth harmonics of this. 
The dashed black line is a fit to the continuum power. The inset shows the weighted photon light curve 
folded on the 16.6-day period. 
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Fig. 2. The orbital modulation of the flux and spectral indices of 1FGL J1018.6-5856 in the 0.1 to 
10 GeV band as measured with the Fermi LAT. r ot and Ti-io are photon spectral indices for 
energies below and above 1 GeV, respectively, using a broken power-law model. 


LS 5039 but is unlike the behavior of LS I +61° 
303. If the high-energy electron distribution 
remains constant along the orbit, spectral changes 
due to the anisotropic IC cross section are ex- 
pected only if the inclination is substantial. In this 
case, harder spectra are expected to occur when 
the stellar photons are forward-scattered by the 
electrons (i.e., at inferior conjunction), which is 
also typically when the scattering rate is at its 
orbital minimum. However, for 1FGL J1018.6— 
5856 the hardness ratio and flux are correlated, 
unlike for LS 5039 (8). If periastron passage co- 
incides with inferior conjunction, then a high pho- 
ton density might compensate for the unfavorable 
interaction angle, but this requires fine-tuning of 
the orbital parameters. The spectral variabili- 
ty is more likely to reflect intrinsic variations — 
for instance, in the cooling of emitting particles. 
Moreover, both PSR B1259-63 and LS I +61° 
303 (7, 9) show that a simple model may not be 
correct The phasing of gamma-ray maximum at 
GeV energies is not consistent with IC scattering 
on stellar photons, as it is delayed in both PSR 
B 1259-63 and LS I +61° 303, implying that 
other mechanisms may be at work. For exam- 
ple, there could be other seed photon sources, 
Doppler boosting, or other radiative mechanisms 
at work. 

The gamma-ray energy flux of 1FGL J1018.6— 
5856 implies a luminosity of ~8 x 10 35 (d ! 5 kpc) 2 
ergs s l (E> 100 MeV), whereas the implied x-ray 
luminosity is highly variable with fluxes up to 
-10 34 (d/5 kpc) 2 ergs s" 1 . For comparison, the 
gamma-ray luminosity of LS 5039 is -2 * 10 35 
(d/2.5 kpc) 2 ergs s' 1 (26). This is somewhat sur- 
prising; relative to LS 5039, the orbital period 
of 1FGL J101 8.6-5856 (longer by a factor of 4) 



Fig. 3. Swift XRT x-ray (left) and UVOT-W1 (right) images of the region around 1FGL J1018.6-5856. The x-ray-optical counterpart is marked by an 
arrow near the center of both images. The LAT 95% confidence ellipses from the 1FGL (12) and 2FGL (29) catalogs are marked. 
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Fig. 4. X-ray (upper paneQ and radio (lower panel) observations of 1FGL J1018.6— 5856 folded on the 
orbital period. The x-ray data are from the Swift XRT and cover the energy range 0.3 to 10 keV. For the 
x-ray observations, the different colors indicate data taken from different 16.58-day orbital cycles. For 
the radio data, green diamonds indicate 9-GHz data and red circles 5.5-GHz data. The radio data are 
from the ATCA. 


implies a major axis that is larger by a factor of 
2.5, so that the mean stellar radiation density 
seen by the compact object is smaller by a factor 
of 6. Tlie higher gamma-ray luminosity of 1FGL 
J1 01 8.6-5856 indicates that the power injected 
in nonthermal particles must therefore be substan- 
tially higher in 1FGL J1 01 8.6-5856 than in LS 
5039. The similarity to LS 5039 suggests that we 
may be observing a rapidly rotating neutron star 
interacting with its companion. This raises the 
possibility that the neutron star rotation period 
might be detectable, as is the case with PSR 
B 1259-63. However, our observations cannot 
definitely exclude an accreting neutron star or 
black hole. 
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Universal Signatures of Fractionalized 
Quantum Critical Points 

Sergei V. Isakov, 1 Roger G. Melko, 2 Matthew B. Hastings 3,4 * 

Ground states of certain materials can support exotic excitations with a charge equal to a fraction of 
the fundamental electron charge. The condensation of these fractionalized particles has been predicted 
to drive unusual quantum phase transitions. Through numerical and theoretical analysis of a physical 
model of interacting lattice bosons, we establish the existence of such an exotic critical point, called XY* 
We measure a highly nonclassical critical exponent r\ = 1.493 and construct a universal scaling function 
of winding number distributions that directly demonstrates the distinct topological sectors of an 
emergent Z 2 gauge field. The universal quantities used to establish this exotic transition can be used to 
detect other fractionalized quantum critical points in future model and material systems. 


I t is noteworthy that, in this age of high-energy 
accelerator experiments, certain types of fun- 
damental quantum particles can only be studied 
in tabletop condensed-matter physics experi- 
ments. Consider, for example, the electron, carrying 
fundamental charge e. Unlike the proton, whose 
charge originates from quarks with fractional 
charge, no energy is sufficiently high to break up 
the charge of an electron. However, as demon- 
strated by the measurement of fractional Hall con- 
ductance, if one places an electron in certain clean 
two-dimensional (2D) materials in a strong mag- 
netic field, its charge can break into fractions — e/3, 
e/5, and so on (/)— with each fractional charge 
arising from a quasipaiticle emerging as an ex- 
citation of the ground state. 

Such quasiparticles share all the important 
characteristics of real particles. In deconfined 
phases of matter with a gap to excitations, the 
quasiparticles can be separated a large distance 
from each other, making them well-defined lo- 
calized objects with a sharp energy-momentum 
dispersion relation. For more than a decade, 
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condensed-matter physicists have searched for 
fractionalized particles in systems other than EM 
effect materials. Theoretical predictions have iden- 
tified a class of low-temperature paramagnets, the 
quantum spin liquids, as holding particular prom- 
ise for supporting them (2- 4). However, exper- 
imental searches for these fractional charges and 
their parent spin-liquid vacuum remain uncon- 
vincing. This is partly a result of the difficulty in 
constructing measurements that are able to iden- 
tify their experimental signatures in the variety 
of materials —ranging from fabricated solid-state 
devices to organic magnets — thought to harbor 
spin-liquid states (5). An interesting recent ex- 
perimental candidate is a set of materials that 
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may display fractional particles with a gapless 
Fermi surface (6-7/), where interaction between 
quasiparticles makes it problematic to even de- 
fine a fractionalized excitation (12). In bosonic 
systems, such interacting gapless fractional quasi- 
particles have been proposed to mediate quantum 
critical points that exist in certain order-to-order 
transitions (13). For these Landau-violating crit- 
ical points (which rely on long-wavelength fluc- 
tuations of fractional particles), signatures of the 
fractionalization are manifest in universal quan- 
tities, such as critical exponents, instead of spe- 
cific material-dependent quantities. 

A major goal of the theoretical community 
has been to demonstrate the existence of these 
“deconfined” quantum critical points in realistic 
microscopic models, with the use of large-scale 
numerical simulation (14). In this work, we used 
quantum Monte Carlo (QMC) to study decon- 
fined quantum criticality in the context of an 
order-to-disorder transition between a superfluid 
and a gapped spin-liquid state, in a physical mod- 
el of lattice bosons. The model we examine is a 
variant of the Bose-Hubbard Hamiltonian, in- 
troduced in (IS), which is a simple model of 
hard-core bosons hopping (with strength t) on 
a 2D kagomc lattice with an energetic constraint 
(V) favoring three bosons per hexagonal lattice 
plaquette:// = ~tj^)[b]bj + bib]] + F£o(«0 f 
(i and j are nearest-neighbor lattice sites; fcjand 
b x are boson creation and annihilation operators. 


Fig. 1. A section of the to- 
roidal kagome lattice sim- 
ulation cell. Gold spheres 
label sites occupied by 
bosons in a representative 
configuration. Red spheres 
are fractional charges, mark- 
ing "defect hexagons" (those 
that do not have three 
bosons per site). Defect hex- 
agons are shown joined by 
an open string. We also show 
a representative dosed string 
not associated with frac- 
tionalized particles. 
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